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We reconsider the well-known and long-debated problem of the calculation 

in 



of the eigenvalues of the Herbst Hamiltonian 2y / p 2 + m 2 — n/r. We give a 
formulation of the problem that allows, for the first time, a perturbative eval- 
uation of the eigenvalues for any n and I, and in principle up to any order in 



k via standard Kato perturbation theory. We present the evaluation of the 
energy of the n = 1 and n = 2 states up to k 6 , confirming the result previ- 
ously obtained by Le Yaouanc et al. with a completely different technique. 
Moreover we give the n = 2, 1 = 1 level, which is new. Discussion of the 
results and comparison with previous findings are given at the end. 
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I. INTRODUCTION 



In the last years the so-called Herbst Hamiltonian (sometimes also called the spinless 
relativistic Coulomb Hamiltonian): 

H = + m 2 - (1) 

has been intensively studied in terms of both the spectrum and the eigenfunctions |T]-§|. 
The reason of the interest in such a Hamiltonian is not only a mathematical one [jTJ but 
mainly a physical one, also due to the relevance of ([I]) in quarkonium phenomenology J7[. In 
fact we recall that, for an appropriate potential V the Hamiltonian H = 2y/p' 2 + m 2 + V(x) 
is the spinless Salpeter equation representing a well-defined standard approximation to the 
Bethe-Salpeter formalism for the description of quark- ant iquark bound states ; in this 

context Eq. (p]) is suitable for the treatment of short-range effects in bound qq systems and 



for the study of decay rates; moreover it was also applied to the study of boson stars [10|. 

For the reasons presented above, a lot of attention has been addressed in the last years 
to the calculation of the ground state E\q of the Herbst Hamiltonian. First, an exact closed 
expression was found in Ref. || which, however, turned out to be wrong (|||; then, many 
analytic estimates of lower p|,[T0| and upper || bounds were given for E w . Up to now we 



are aware of only a paper devoted to the perturbative calculation of E w ||. The method 
used in || has the merit of being systematic and allows the authors to settle an iteration 
procedure on a reasonably rigorous basis Q, but it is not general. Its starting point is a 
peculiar representation of the Hilbert space of I = states (introduced in 0) in which the 
Coulomb potential operator has a particularly simple form. 

For this reason we have decided that it could be useful and instructive to apply a general 



1 The authors of Ref. Q derive an equation for a function related to the eigenfunction and use it 
to obtain contributions up to k 5 to the I = levels. Then, they use the analyticity properties of 
the wave function to obtain next to leading orders up to k 7 . 
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and powerful perturbative technique to the calculation of the eigenvalues of the Herbst 
Hamiltonian. This method, not being founded on peculiar transformations, can in principle 
be applied for any value of the angular momentum and also extended to the consideration 
of potentials different from the Coulombian one. In fact it was proposed for the first time 
for the calculation of the positronium energy levels in Ref. 

The plan of the paper is as follows. In Sec. II we explain the perturbative technique; 
in Sec. Ill we apply it to the evaluation of the energy levels in the case n = 1 and n = 2 
(7 = 0,1); in Sec. IV we make a comparison with the results previously obtained in the 
literature and draw some conclusions. 

II. DESCRIPTION OF THE METHOD 

We have to solve the equation 

= E<& (2) 

with the two-body relativistic Coulombian Hamiltonian given in (|I|) or equivalently with the 
one-body Hamiltonian 

H = ^ p 2 + m 2_ ^. ( 3 ) 

The energy levels of the two-body (]l|) and the one-body @ cases may easily be related by 
identifying both mass and Coulomb strength coupling parameters according to m — > 2m 
and a — > k/2. In the following we consider the one-body case. 

The perturbative solution of Eq. (0) which comes immediately to mind refers to the 
well-known non-relativistic Schrodinger Hamiltonian as the unperturbed starting point and 
supplies the eigenvalues E nt (n = 1, 2 . . . ; I = 0, 1, . . . , n — 1) by evaluating perturbatively 
the expectation value of the relativistic corrections, coming from the expansion of the kinetic 
square root, on the analytically available zeroth-order eigenfunctions. This can be simply 
done for the contribution in the fourth power of p. For instance, for the ground state 
zeroth-order energy E w , we can calculate 
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5E 10 = (ip w0 \ - ^-t|v?ioo) 
Sm 6 



64 



-ma 



(4) 



where <^i o is the ground state Schrodinger-Coulomb wave-function: 

8^(ma) 2 



^iooIPJ 



IP' 



ma 



(5) 



However, already with the sixth power of p, one encounters a divergent integral; in order 
to get the next term in the perturbative expansion, which turns out to be of order a 5 , one 
has to sum up an infinite number of contributions coming from all the possible intermediate 
states (cf. Ref. [[]]). In general, the relativistic kinetic energy a/p 2 + m 2 , behaving as |p| 
for very large |p|, cannot be approximated by any finite polynomial in p 2 . This seems to 
prevent the possibility of profiting from the known results for the Schrodinger equation in 
the relativistic case. There is, however, a trick for recovering, at least formally, a "kinetic 
energy" p 2 , even starting with ^P 2 + m' 2 ; the price to pay is a kind of energy- dependent 
potential instead of the Coulomb one [11]. Let us consider Eq. (|]) written in the momentum 
space 

d 3 q 



(2tt) 3 
V(P) = 



(2tt) 3 5 3 ( P - q)( y/p 2 + m 2 - E) + V{p - q) 

4na 



$(q) = 0. 



(6) 



We define 



2m 



E v + E 



E,, 



p2 _|_ m 2 



(7) 



and put 



$(p)=i?- 1 ( J E;;p)0(p) 
V(p - q) = R(E; p)v(E; p, q)R(E; q). 

By using Eqs. (^) and (||), we rewrite Eq. (|5]) in the form 



(8) 



At least, not without some ad hoc transformations and intricate calculations as in Ref. H 
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d 3 q 



(2tt) V(p - q) 



£ 2 -m 2 



v(£;p,q) 



(q) = 



(9) 



2m 2m 

which is formally of the Schrodinger type, i.e. with a quadratic kinetic term, but with 
an energy-dependent interaction. Equation (||) is the appropriate starting point for the 
application of a perturbative technique. 

Let us rephrase the same transformation in the propagator formalism, which is more 
appropriate to treat an energy-dependent perturbation term. With reference to Eq. (0) we 
can write for the one-particle propagator G : 



G(£;p,q) = G (£;p) 
Go being the free propagator 



r d 3 k 

(2n) 3 5 3 (p - q) + J j^V(p - k)G(E; k, q) 



(10) 



E — E„ 



(11) 



Now, by using the identity 

1 



1 



( 



E + E p 

E-E p (E 2 - m 2 ) /2m - p 2 /2m V 2m 
Eqs. (H), and defining 

1 E 2 — m 2 



(12) 



9o(E* 



E* - p 2 /2m 



E 



2m 



(13) 



we can rewrite Eq. (|T0|) as 

g(E*;p,q)=g (E*;p) 
where v is given in (||) and 



r H 3 k 

(27r) 3 <5 3 (p-q) + J — t;(£;p,qM£*;k,q) 



(14) 



g(E*; p, q) = R{E; p)G{E; p, q)R(E; q). 



(15) 



We emphasize that in Eq. ( |T^ ) the propagators g and go has to be understood as depending 
on the modified energy E*, which is a function of E. 

Inserting V in place of v in ( |T1D we obtain the exactly solvable equation for the well- 
known Schrodinger-Coulomb propagator g s : 



9s(E*;p,q) = go(E*;p) 



r r/ 3 k 

(27r) 3 5 3 (p - q) + J j^r- 3 V(p - k)g s (E*; k, q) 



rf 3 k 

The g s propagator has poles in correspondence of the energy values 
ma 2 „ „„ i — ma 



(16) 



E* = -— ^E = E» = ^m 2 -rt, = — , n = l,2,... , (17) 

and can be written for any n as 

; P, q = ^— J77T1 — + 9n{E ; p, q), 

hfi- n (^ ;p,q) 



m E/,m^nim(p)<^/ TO (q 



K E -K 



+ g' n (E;p,q), (18) 



V'nim being the Schrodinger-Coulomb wave-functions and g n (or g' n ) being the regular part 
of the propagator. 

Now, by using standard Kato perturbation theory fl2|, one obtains the following expan- 
sion for the energy levels (5V = v — V), 

E nl = E° n +^- [ (5V(E n )) nl + 

+ ^{mKmK)mK))ni + ^(m^u (^^m) } + 

+ ^((5VW B M^(W B (^)^(^)>^ + ... , (19) 
where the symbol ( ) n i stands for 

i)nl = K^Jj2^Jj2^( )^ m (p)^ m (q). (20) 
The summation over the quantic number m and the factors 1/(2/ + 1) takes in account the 
degeneration of the level. 

Expansion ( |T9| ) is valid for any n and any I. As higher order in 5V supply higher-order 
leading contributions in a, Eq. ([19]) is perturbative not only in 5V but also in a. So up to 
a given order in a only a finite number of terms in ( |I9"D contribute to it. 

In the next section we apply Eq. ( JT9] ) to the evaluation of E n i in the cases n = 1 and 
n = 2 I = 0, 1 up to the order a 6 . 
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III. CALCULATION OF THE ENERGY LEVELS 



Up to order a 6 , only the terms explicitly written in Eq. ( M ) have to be taken in to 
account. It is useful to write g'(E®) in the form: 



2 ' 



2m J 2E* 1 i 



nlm 



rn 



where 



4ifl + *(^)Mif) +J ^ (21) 



K( E n) = R n - 777=^2 VnlmVllm- 



2m J 2Er f 

The functions R\ (— ji/2m) and i?2 ( — 7l/2^) are given in the appendix. From Eq. fl2~T|) 
and 

(*v(B»)9;(B»)^(B»)j;(£;2)^(E, , ;))„,= 

we obtain the only relevant contributions to the energy levels up to order a 6 : 

E "< = E ° + 2FTT™>- + 27Tt{ (if) ^ B °>)„, + 

+ ^V(£j>,) S „ Vjo iV(B„ )^ + (5V(K)R'n(K)6V(E „))nl + 

+ arr(^* (if) ^ (£ ° )9 » (if) ^ (£ ° ) )„ 1 +o(a8) - (22) 

With the symbol o(a; 6 ) we indicate higher-order contributions, typically starting with a 7 
and a 7 In a terms. 

First let us consider the case / = 0. Using the Schrodinger equation (goVip = if) we can 
perform some straightforward cancellations and reduce Eq. (p2|) to the form: 

+ mK)K(E° n W(E Q n )) n0 + o(a 6 ). (23) 
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For the level n = 1 we obtain: 



{vg vg Q 6V) 1Q 1 4 8 a 5 6 , 

= — a H a ma + 

m 2 3 7r 



Hc(3)-^+>>-^w + °< a6 >- 



m 4 



Inserting the various contributions in Eq. (p3|) , we obtain the final result: 

El ° 1 - -a 2 - -a 4 + |— - a 6 In a" 1 + f ^(3) - ^ - - V + «(« 6 )- (24) 



to 2 8 3 vr Vtt 2 ^ 7 7T 2 16/ 

Similarly for the n = 2, 1 = case we calculate 
(^0^0(^)20 3 4 la 5 1 6 



TO 

11 7 1 w „, 7T 2 I 

4t^ + 8^ 2 " ,v "" / 32 8 



+ - oC(3) - 7 ^ + -^C(3) - — - - In 2 U 6 + o(a 6 ), 



(SV) 20 (^5V) W 



m 

(SVR' 2 5V) 20 
m 



a" 

32 


+ 3~- 


/127 
\512 




—a 8 
256 


+ o(« 6 ), 


7T 

V32 


231 


~ 512 + 



+ ^C(3))a 6 + o(a 6 ), 



and then obtain 



£20 , 1 2 13 4 1« 5 1 6 , 

= 1 a a H a ma + 

to 8 128 3 vr 8 

+ G>»-ii-i^-i ln2 ) a8+ ^' < 25 > 

When Z is different from zero, we have to apply directly Eq. (0). Let us consider the 
case n — 2, 1 — 1. Then 

(5V) 21 1 4 19 6 . 6 , 



to 32 512 



(^0^)21 / 1 2 15 

= — 7T a + o(a 

in V96 128/ V 

— = — vr H — C 3 cr + o(a ), 

V48 8 SV ; 256/ V ; ' 



/// 



(^^) 21 = /2051 _l ± A 

m V4608 8 SV ; 32 / v ; 
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(SV) 21 (^6V) 



121 



a 6 + o(a u 



m 



256 



(6Vg SVg SV) 21 



m 



o(a 6 ). 



and, summing up all the contributions, we obtain 



E21 1 2 7 4 25 

= 1 a ol H 

m 8 384 27648 



a 6 + o(a 6 ) 



(26) 



which, surprisingly, does not contain contributions in a 5 and a 6 In a. 

The description (in another context) of some techniques used to extract from the integrals 
the various order in powers of a can be found in ||13|| . All calculations were performed with 



the help of the program of symbolic manipulations FORM [14 . 



IV. DISCUSSION AND CONCLUSIONS 



In conclusion, in the case n — 1, 2, I — 0, we have obtained a result (Eqs. (24) and 
(P5|)) that confirms, up to order a 6 , the previous findings of Le Yaouanc et al. |J. However, 
the result obtained in the case n = 2, 1 = 1 (Eq. (^)) is new and displays also a different 
analytical structure. 

We emphasize again that our method is general and can be applied to the calculation 
of the eigenvalue of the Herbst Hamiltonian for any n and any /, and in principle up to the 
desired order in a. Moreover, due to its generality, this method could be applied also in 
the case of a potential different from the Coulomb one. We notice that the problem can be 
solved on a general ground and by well-known energy- dependent perturbation theory due 
to the appropriate formulation of the starting point in ([|). In this way, in fact, we gain an 
energy-dependent perturbation containing a. 

Finally let us mention some other points: 

1) The result we have obtained for the ground-state energy fits nicely in the upper and 
lower variational bounds given in Refs. p,IOfl. 
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2) The perturbative expansion of the eigenvalues given in (p^)-(pq) is suitable to obtain 
good numerical estimates of the energy of the level in principle up to the critical value 
of the coupling constant a = 2/n. As an example, we present in Table 1 the evaluation 
of E w , E 2 o and £21 for different values of a. We have chosen the same values of a as in 
Ref. [ID], so that it could be appreciated (for small couplings at least) the improvement 
introduced by this perturbative expansion with respect to the variational evaluation 



(cf. the variational upper E™ x and lower E™ n estimates given in Table 2 of Ref. [10 



and transcribed here in columns 2 and 3 of Table 1). In Table 1, the last figure in 
our results has been put in parenthesis because it is affected by the error coming from 
the neglected a 7 and a 7 In a contributions. For a > 0.5 the error obviously grows 
and we have therefore not considered significant to present the results and to make a 
comparison also in this case. 



Acknowledgements 

We would like to thank Wolfgang Lucha for useful discussion and encouragement, and Ste- 
fano Laporta for the evaluation of some integrals. We are extremely grateful to Andre 
Martin for pointing to our attention an error in the final result for the n = 2, / = 1 state. 



APPENDIX A: 

The non-singular part g n of the Schrodinger-Coulomb propagator is obtained by expand- 
ing the analytical expression of g s in the vicinity of the energy level and by subtracting the 
pole contribution. The expression of g\ evaluated in —7^ /2m reads (cf. ||T5"H ): 

~ /-7i 2 \ (2tt) 3 <5 3 ( P - q) 1 , 1 



2m' ' J — 7i/2m — p 2 /2m — 7i/2m — p 2 /2m — 7 2 /2m — q 2 /2m 

647rm7 3 (\ _ { 7i _ 4 _ ll \ 



(p 2 + 7 2 ) 2 (q 2 + 7 2 ) 2 \2 p 2 + 7l 2 q 2 + 7l 2 
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64^7717? f 1 2Ci - 1 



lnCi + ' arctg v /4Ci - 1 



(p 2 + 7?) 2 (q 2 + 7i 2 ) 2 I 2 v/iCW 

In a similar way one obtains the regular part corresponding to the n = 2 pole, which 
evaluated in — 7 2 /2m reads: 

/- 72 2 \ (27r)^ 3 (p-q) 1 1 
#2 ; P, q = 277^ 777^— + 2 77^ T77T— ^(P ~ 0) 



2m' ' y — 7 2 /2m — p 2 /2m — 7|/2m — p 2 /2m — 7|/2to— q 2 /2m 

2567rm 72 3 [ 2 7 |(p - q) 2 / 9 , 6 7 | 6 7 | 



^P 2 + 7l) 2 (q 2 + 7 2 2 ) 2 I (P 2 + 7 2 2 )(q 2 + 7 2 2 ) V" 2 + p 2 + 7 2 2 + q 2 + 7 2 2 J + 



+ 



3 4 7 2 4 72 2 



2 p 2 + 72 2 q 2 + 72 2 
( 1 ^ , 2C 2 -4+l/C 2 



with 



r (P 2 + 7 2 )(q 2 + 7 2 ) 



' II 



47n(p - q) 

and 7„ given in Eq. (|T^). 
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TABLES 



TABLE I. Calculation of the energy levels given in Eqs. (24), (25), (26). In columns 2 and 3 
are shown the minimum E™ 1 and the maximum E^ x variational estimates of the ground-state 
energy as given in Table 2 of [10] are shown. The figure in parenthesis is of the same order as the 
uncalculated a 7 and a 7 In a contributions. 



a 


Ef$ n /m{W] 


E^ x /m[W] 


E w /m 


E 2 o/m 


E 2 i/m 


0.0155522 


0.9998785 


0.9998791 


0.99987902866(7) 


0.99996976027(9) 


0.99996976506(8) 


0.1425460 


0.989458 


0.989613 


0.98960(4) 


0.99742(0) 


0.99745(2) 


0.2599358 


0.96309 


0.96364 


0.9635(1) 


0.9911(0) 


0.9914(7) 


0.3566678 


0.92578 


0.92673 


0.926(1) 


0.982(3) 


0.983(8) 


0.4359255 


0.88013 


0.88139 


0.88(0) 


0.97(2) 


0.97(5) 


0.5000000 


0.82758 


0.82910 


0.82(9) 


0.96(1) 


0.96(7) 
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